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Abstract We study p — w— interference by analyzing J/ip — ► 7r + 7r~7r 7r . PDG-2002 data 
on J/tp decays into PP and PV {P denotes pseudoscalar mesons; V, vector mesons) 
' are used to fit a generic model which describes the J/tp decays. From the fits, we ob- 

t^J" ■ tain anomalously large branching ratio Br(p° — * ir + Tr~ir°) ~ 10~ 3 — 10~ 2 . A theoret- 

C> ical analysis for it is also provided, and the prediction is in good agreement with the 

anomalously large Br(p° — > 7r + 7r _ 7r°). By the fit, we also get the rj — ry'— mixing angle 
O 6 = —19.68° ± 1.49° and the constituent quark mass ratio m u /m s ~ 0.6 which are all 

q^. reasonable. 

<D '. 

' Keyword J/ip— decays, p — lo interference, SU(3)-breaking effect, 77 — 77' mixing 

>: 

^ ! 1 introduction 

.5?. 

Recently it has been predicted in refill that there is a large isospin symmetry breaking 
enhancement effect in the decay p° — > 7r°7 comparing with — > 71^7 due to p — w-interference, 
which was called as hidden isospin-breaking effects in HI. This prediction has been confirmed 
by the renewed data in PDG-2002^. Following the discussion of LU, it could be expected 
that a similar large hidden isospin-breaking effect should also exist in p° — ► 3tt due to p — u> 
interference. In this paper we try to analyze J /ip — > (p , u))(tt°, 77, 77') — > (37r)(7r°, 77, 77'), and to 
reveal p — w-interference effect and then finally to abstract out the branching ratio of p° — ► 3tt. 
Such an analysis should be necessary for further confirming the enhancement effects mentioned 
above, and be also interesting for the G-parity violating process studies. 

It is very difficult to directly measure B r (p° — + 37r) experimentally both because T p >> 
(rria, — mp) and because the G-parity conserving decay mode of uj — > Stt is dominate. This is 
the reason why there is still no a reliable value for B r (p — > 37r) yet so far in the literatureElEl , 
Fortunately, this quantity can be obtained by fitting the data of J/ip — > PP and PV (where P 
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denotes the pseudo-scalar meson nonet (tt, K,rj,r]'), and V, the vector mesons (p, uj, K*, (p)) <z*. 
Actually, more than fifteen years ago, B r (p° — - > 37r) was estimated in ref.El by using the MARK- 
III data of (J/ip — > PP, PF)-decays with P = (ir,K) and V = (p,co,K*). However, people 
including the authors of ref.El does not think their result of B r (p° — > 37r) is very reliable (see the 
discussion in ref.Sl and ref.El). The reasons are multi-ply, and some of them may be as follows: 
1, the J j ' ■(/'-decay data quality at that time was not good enough ; 2, the fitting is not complete 
because the processes of {J/ip — > rpV, rj'V, P<j>) were not considered; 3, lacking a theoretical 
understanding why their result is so significantly different from the result of ref.^ which is 
quoted by pdgEI. Our motive of this paper is to try to solve those problems: 1, We shall use 
nowadays data of ( J/ip PP, PV) in PDG-2002EI to perform the fit; 2, (J/ip ->■ 77V, rfV, Pep) 
will be considered in our analysis; 3, And the rationality of the result will be argued, i.e., we'll 
see that the result is just consistent with the theoretical analysis in ref.d. 

The contents of the paper are organized as follows. In section 2 we describe the p — 
uj— interference in the process J/ip — > 7r + 7r _ 7r°7r°, and give the the branching ratio formulas 
for (Jip — > PP, PV). In section 3, by using the PDG-2002 J/ip decay data and the formulas 
given in the section 2 we perform the datum fits. The Br(p° — > ir + Tr~ir°) is obtained. The 
section 4 is devoted to estimate Br(p° — > 7r + 7r~7r ) by theoretical analysis. Finally, we briefly 
discuss the results. 



2 p — uj Interference and branching ratios of J/i/j decays 
into PV and PP 



The elusive p° — > tt + tt 7r° decay could be observed in J/ip decays into the tt + it 7r°7r° 
final stated. Indeed, this decay can proceed through the interfering channels J/ip — * luit — > 
7T + 7r~7r 7r° and J/ip — > pir° — * 7r + 7r~7r 7r . Because J/ip — ► pir° is caused both by strong 
interaction via 3 gluons and by electromagnetic (EM) interaction and J/ip —> luit is caused 
by EM interaction merely, the passibility of the decay J/ip — > pir° is much larger than one of 
J/ip -> uir°, i.e., Y(J/ip -> pn°) » T(J/ip -► ujtt°) (by using the 2002-PDG data^l we have 
T(J/ip — > pn ) ~ 2.5 x 10 2 T(J/ip — * W7T )). Consequently, even though T(oj — ► 7r + 7r~7r°) may 
be much larger than T(p — * 7r + 7r~7r°), it is still hopeful to measure T(p — > 7r + 7r~7r°) by studying 
the p — w-interference effects in the process of J/ip — > 7r + 7r~7r 7r . 

To J/ip — > (p,uj)ir° — ► 37T7T , the corresponding s-dependence Breit-Wigner is written as^ 
F(s) = BW u {s) + ee l9 ' BWp{s), where 



W = 2 1 ■ r (1) 

is the normalized Breit-Wigner curve for i = lu, p resonance with the mass of mi and total 
width of r,*. The factor e(O') is the modulus (phase) of the amplitude proceeding through the p 
resonance relative to uj resonance, which can be written as: 



L4(J/V> ^ A°)l ^ lBR(p0^3n) ^ 



\A(J /ip — > a;7r°)| y BR(oj -> 3tt) 
and 0' = 0^/^-,^) + %^3tt) witn * = P, UJ - 
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As we have the relations m p ~ m w = m and T p =/= T^, the total effect is integrated over s 



ds\F(s)\ 2 = l + e 2 + 2ecos#'^^A (3) 



2 /r r 

The third term of the expression 2e cosfl'( p- P ), is the p—u> interference term. According 

to ref.Bl 1 , the phase 9' is equal to zero, and the interference effect produces a magnificent factor 
as a whole(we assume that the w resonance contribution is one),ie. 



ds\F{s)\ z = l + e 2 + e, (4) 



where T p ~ 16r w l^l has been used. Then, the interference between p and uj in the J/tp — > 47r 
provides a relation as follows 

BR(J/tp -> w(p)tt° -> 4vr) = (1 + e 2 + e)BR(J/ip -► wvr° -> 4tt). (5) 

The BR(J/tp — ► uj(p)iv — > 47r) can be detected directly in experiments, but there is no a direct 
experiment way to measure BR(J/ip — > low — > An) in the right hand side of the above relation. 
Fortunately, it can be got by fitting the branch ratios of J/tp — > PP and PV( where P and 
denote pseudoscalar- and vector mesons respectively) El. As both BR(J/tp — > Lo{p)ir® — > 47r) 
and BR(J/tp — > oj7r° — » 47r) are known, we will have e by Eq.JSJand then obtain desired quantity 
PP(p° -> 3tt) via Eq.HJ 

Decays of J/^ into (PV) and into (PP) can be factorized by a very simple and general 
consideration described as follows^. The decays proceed through a strongly interacting three- 
gluon (ggg) intermediate state and through electromagnetic interaction mediated by one photon 
7 and (jgg) states. In the gluonic case, there are two 1 = transitions cc — > ggg — ► LHH+^S an( j 
cc — > — > ss, which are proportional to the amplitude A and ^= respectively. The parameter 
A accounts for flavor-iSC/ (3) breaking effect. The electromagnetic transitions generate the 1=1 
state ( uu ^ d ) a nd two 1 = states ( ( uu ^ d ) anc [ S g), Their amplitude are proportional to 3a, a 
and — ^/2\a respectively. The flavor space wave functions for P including 77 and 7/ read 



p = A a $ a = v/2 



(6) 



where 77s = rj cos 9 + 77' sinf? and 771 = 77' cos# — 77 sin 9 with 9 as (77 — ?7')-mixing angle, and the 
1^-wave functions including </>'s read 

p+ K*+ 

V v2| -^_+^ iY*° I- IT) 

K*° 

With the above, the decay amplitudes of {J/tp —> PP, PV) are as follows 




A(tt + tt-) = 3a, (8) 



1 9 = in J/tp —* (p or a;)7r° — > 7r+7r 27r° and 8 = tt/2 in J/ip — » (p or w)7r° — > 7r+7r 7r° in the paper of 
A.Bramon and J.Casulleras,Phys.Lett.,1986,173B: 97. 
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A(K+K-) = ^(l-X)A + (2 + X)a, (9) 

A(K°K°) = ±(l-\)A-(l-\)a, (10) 

A(p\°)=f v (A + a), (11) 

A{K*+K-) = + \)A + (2 - X)a], (12) 

A(^ /?«) = f v [\{l + X)A - (1 + A)a], (13) 

A( W7 r°) = f v (3a), (14) 

A(p»/) - fv(3aX v ,), (15) 

A(wr/) = /„[(i4 + a)X n , + V2rA(V2X v , + Y v >], (16) 

Ml«l) = fv(?aX v ), (17) 

A(urj) = f v [(A + a)X v + V2rA(^2)X v + Y v ], (18) 

A(<f> V ) = f v [(A - 2a)XY v + rA(V2X n + y,)ii±% (19) 

MH) = fv[(A - 2a)\Y v , + rA(V2X n , + fy) (I±*l], (20) 

where X v — ^/i cos# — y|sin0, X v i = ^Jlj sin 9 + «^|cos#, Y v = —X v i,Y v i = X v . The 
additional parameter r is the relative weight of the disconnected diagram to connected diagram 
for the decays involving the final state rj or The Eqs.|SJ-^]are same as ones in ref.^, 

and others are new. 

The corresponding branching ratios of these decays are following 

Br{n + TT-) = 9a 2 , (21) 

Br{K+K-) = \hl - \)A + (2 + A)ae^| 2 , (22) 

Br(K°K°) = |i(l - X)A - (1 - A)ae'^| 2 , (23) 

Br(p°Tr°) = f v \(A + ae^)\ 2 , (24) 

Br(K*+K-) = / 2 |i(l + X)A + (2 - A)ae^| 2 , (25) 

Br(K*°K°) = / 2 |i(l + X)A - (1 + A)ae^| 2 , (26) 

J Br(c(p ) 7 r°) = (l + e + e 2 )/> 2 , (27) 

Br(pr/) = f*\3aX v ,\ 2 , (28) 

Br(Lu V ') = f 2 \(A + ae^)X v , + V2rA(V2X v> + Y v ,\ 2 , (29) 

Br(p V ) = / 2 |3aX„| 2 , (30) 

Br{Lo V ) = f 2 \(A + ae^)X v + V2rA(^~(2)X v + K„| 2 , (31) 
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Br{H) = %\{A - 2ae^)XY n + rA(V2X v + Y n )—^-\ 2 , (32) 

Br(H') = fv\(A - 2ae^)\Y v , + rA(V2X v , + Y v ,)^-^\ 2 , (33) 

where <fi is their relative phase between A and a, and the parameter A and a are real. In the 
Eq. |31 the p — uj interference effect is subtracted from the branching ratio of J ftp — » ptt°. In the 
Eq. 123 a magnificent factor 1 + e + e 2 has been added due to Eq. [5] in order to taking the p — u) 
interference effects into account. Actually, through directly detecting the data of J/tp — > An 
one can only get Br(w(p°)ir°) rather than Br(um°) which is equal to / 2 9a 2 . 

In the branching ratio formulae of Eas. l21l - I33l we do not write out the corresponding phase- 
space factors explicitly which are proportional to the cube of the final momenta in two-body 
decays. They will be taken into account in the practical phenomenological fit later. 



3 Datum fit to obtain BR(p — > tt + tt tt°) 

Let's now use the data in PDG-20020 to perform the fit to all branching ratios of Ea. l21ll3*3l 
This fit will lead to determining the e and BR(p — > 7r + 7r~7r°). The experimental branching 
ratio data of PDG-2002 are listed in the second column of Table 1. 

Firstly, following ref.EI, we use the branch ratio data of (J/tp — ► PP, V(tt°,K)) only to 
perform a fit to Eos. 1211- I27f call it as (PP, V(tt°, K))— fit hereafter). In this case, there are 
seven equations with six adjustable free parameters a, A, A, <f>, fy, e, and hence it is an 
over-determination problem with potential of predictions. The fit with minimum \ 2 = 0.46 
leads to the values of the parameters and seven corresponding branching ratios listed in the 
third column of Table 1, in which (w7T°) uncor and (unr°) cor represent Br(J/ip — > w(p)7r° — > Air) 
and Br(J/ip — > lutt — > Air) respectively, i.e., 

Br(J/4> ^ un°) cor = / 2 9a 2 . (34) 

In the fit (see Table 1), we have (w7r°)„„ cor = (4.2±0.61) x 1CT 4 , a = 0.21±0.02, A = 2.94±0.72 
and the interference factor e = 0.71 ± 0.58, then we obtain (w7r°) co ,. = / 2 9a 2 = (1.89 ± 0.83) x 
10~ 4 . In other hand, from Eqs.QdHHI), e reads 



\A + ae*\ Brjp^to) 
€ ~ \3a\ X \J Br (cj->3n)- (6b) 

Then the branching ratio of p — > ir + ir~Tr° is predicted as follows 

Br(p -> 3vr) = Br(u - 3.) x ( (36) 

Substituting the a—, A— and 0— values obtained from the fit (see the third column of Ta- 
ble 1) and experiment data of Br(ui — > 37r) into Eq. 02] and Eq. |33 we then obtain the 
(PP, y(7T°, iT))-nt's results as follows 

Br(J/i,^Lu^) cor \ {PPy ^^ K)) = (1.89 ± 0.83) x 10- 4 , (37) 
Br(p^7T + TT-7T )\ iP p iV(w0iK)) = (2.0 ±1.64) x 10~ 2 . (38) 



5 



TABLE 1: The second column displayed the experimental values for the branch ratios 
of J/ip -> PP and J/tp -> PV in PDG-2000 datum . The results of a fit to the 
first seven branching ratios is listed in the third column. The results of a fit to the 
total thirteen branching ratios is listed in the forth column. 



J fyjaecay 


PDG-2002(10 ) 


a partial fit 1(10 ) 


1 1_ 1 p j o / -i r\ — 4 \ 

a global nt2(10 ) 


1.7T 7T 


1.47 ± 0.23 


1 A A 1 C\ C\ O 

1.44 ± 0.23 


1.92 ± 0.05 


2.K + K~ 


2.37 ± 0.31 


2.45 ± 0.28 


2.04 ± 0.08 


3.K K° 


1.08 ± 0.14 


1.06 ± 0.14 


0.87 ± 0.05 


A 

i.p 7T 


42.0 ± 5 


43.04 ± 4.48 


41.97 ± 0.68 


5.K* + K~ 


25.0 ± 2.0 


24.11 ± 1.41 


23.64 ± 0.5 


6.K K 


21.0 ± 2.0 


21.66 ± 1.7 


24.21 ± 0.49 


7 (, 

I . yuJTi )uncor 


4.2 ± U.O 


4.2 ± 0.61 


4.2 ± 0.2 


8 ( 

°-{PV ) 


l.UO it U.lo 




n 7 4- n op; 


9.(aV) 


1.67 ±0.25 




1.73 ±0.11 


10.(pij) 


1.93 ±0.23 




1.82 ±0.08 


U.(urf) 


15.8 ±1.6 




18.32 ± 0.36 


12.(077) 


6.5 ±0.7 




5.85 ±0.23 


13.(0??') 


3.3 ±0.4 




2.55 ±0.23 


x 2 




0.46/1 


21.4/5 


EDM 




0.45£ - 06 


0.69£ - 06 


fit a 




0.21 ±0.02 


0.24 ±0.012 


A 




2.94 ±0.72 


2.69 ±0.17 


A 




0.6 ±0.1 


0.62 ± 0.03 


4> 




1.37 ±0.14 


1.6 ±0.11 


fv 




1.26 ±0.36 


1.38 ±0.1 


e 




0.71 ±0.58 


0.3 ±0.16 


6> 






-0.343 ± 0.026 


r 






-0.144 ±0.001 






(1.89 ±0.83) x 10~ 4 


(3.02 ±0.2) x 10" 4 


p — ► 37T 




(2.0 ± 1.64) x 10~ 2 


(0.59 ±0.315) x 10~ 2 



Secondly, we perform more complete datum fit in which the processes of J /ip — > Vr] and 
J/ip — > Vrj' are included. In this case, there are 13 equations H21I33|1 and eight free parameters: 
a, A, A, (f>, fv, e, 0, r. And hence it is an over-determination problem with more constraints, 
and will be called as (PP, V(ir°, K, 77, f]'))— fit hereafter. The results are as follows 

Br{J/^^LJTT°) cor \ { pp iV{ ^ >KiTliril)) = (3.02 ± 0.2) x 1(T 4 , (39) 
Br(p^ 7r+7r-7r )| (P p : v( T o^^ ir)0) = (0.59 ± 0.315) x 10" 2 , (40) 
6 = -0.343 ± 0.026 = -19.68° ± 1.49°, (41) 

where 77— 77'— mixing angle 9 is agreement with one in ref.ElEl, and both Br(J/i> -> W7r ) cor and 
Br(p — > 7r + 7T _ 7r°) are reasonable agreement with the results l|37l) l|38(l obtained by (PP, V(ir°, K)) — fit 
within the errors. 

The parameter A is the constituent quark mass ratio m u /m s which should be about 0.6^^^ 
due to light flavor SU(3)-breaking. The results of A ~ 0.6 ± 0.1 for (PP, V(ir°, K))-&t and 
A ~ 0.62 ± 0.03 for (PP, V(tt°,K, 77, 77'))— fit indicate the fits meet this requirement, and, hence, 
the results yielded by them are rather reliable. 
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4 Large isospin breaking effect in decay p° — > tt + tt tt° 

In this section, following ref.^, we provide a theoretical estimation to Br(p° — > 7r + 7r~7r°). 
Using Feynman propagators method, the on-shell amplitude^ of the decay p — ^ 7r + 7r~7r° is 
determined by 

M p o^ = (/ p3 . + 2 Upu{ ?l U \ ) 2 , (42) 

where, the momentum-dependent p° — uj interference amplitude Tl puJ (q 2 ) is defined by the p — lu 
interaction Lagrangian £ puJ as follows 



C pu = / * e - i *-*IL pu tf)(fr ^-)^(p)p°Jx). (43) 



(27r) 4 p £ 

The first and second term of expression (|42|) correspond to the contributions of direct coupling 
(p — 3-7r) and w resonance exchange respectively. Because m p ~ m w and r w is small, the 
denominator of the second term is small. Therefore, contribution from u> resonance exchange 
is large. This is called "hidden isospin symmetry breaking effect" according to This effect 
brings a significant contribution and plays an essential role in the decay p — > 7r + 7r~7r°. So when 
we deal with the decay of p — -> 7r + 7r~7r°, the process p — > uj — > 7r + 7r~7r° must be considered. 

In fact, the contributions of p° — uj interference are dominant and the direct coupling can 
be omitted. The direct coupling f p ^ oc (m<j — m u ), therefore, it is very small. In order to 
be sure of this point, we derive this quantity in a practical model called as U(2)l x U{2)r 
chiral theory of mesons ^ in follows. Denoting the direct vertices of p° — 3ir as £ P 3ir — 
J ? p37i-e MiyQ,S 'tijkP^duir' d a ir J dpir k , then / p37r can be calculated in this theory^ and has the form 

rrid — m„ / 16c 6c 2 8c 3 \ „ , , 

^ = -^iiii} ~^ + v-^y- 2x icriiMev - 3 < ™ 

where the values of model's parameters m, g, c determined in ref. SI have been used. To the 
second term in the parentheses of expression l)42fl . H pu (m 2 ) has be determined to approximate 
—4 x 10 3 MeV 2 EHEl. (jj _ > 3^ i s the dominant channel for uj— decays, and hence can be 
estimated by using the width T u ^ n — 7.5MeV. It's approximate value is about 3 x 10~ 7 MeV~ 3 . 
Thus the typical value of the second term in expression l|42l) is (5 + 2i) x 10~ 8 MeV -3 approxi- 
mately. Comparing it with expression l|44l) . we can see that the direct coupling j p ^ is indeed 
very small, and it is ignorablc. Therefore, discarding j ' p ^ in expression 142|) . we have, approx- 
imately, 



n P ^(TOp) 



to 2 — m 2 + im u T u 



2 



(45) 



This equation means that the contributions due to p — a;— interference to Br(p° — > 37r) are 
dominate, or the hidden isospin-breaking effects introduced in ^ are dominate for the process 
p° — > 37r. From expression i|45|) we obtain desired result as follows 

BR(p° -> 3vr) ~ 0.2 x 1CT 2 . (46) 

Our experiment datum fitting result 14U|) is consistent with this theoretical estimation result. 
This fact indicates that both (PP, V(tt°, K))-&t and (PP, V(ir°, K, rf^'^-fit are reasonable 
even though the resulti ng B R,(p° — > 37r) is much larger than one in ref.l^l and rather closes the 
upper limit for it in refM^*. 



7 



5 Discussion 



Through the study presented in the above, we conclude that p — to— interference effects can be 
detected in the J/ip — > ir + ir~ir Q ir° decay, which receives a contribution from the p° — ► 7r + 7r~7r° 
decay mode. J/ip decays offer an almost unique opportunity for observing p° — > 7r + 7r~7r°, 
where the smallness of Br(p° — > 7r + 7r~7r°)/i?r(aj — > 7r + 7r~7r°) is compensated by the large 
ratio A(J/tp — ■> pir°)/A(J/ip — > oj7r°) between a (simply Zweig-forbidden) strong amplitude 
over an EM one. This is the key point for the practical determining Br(p° — > 7r + 7r~7r°) 
through employing J ftp decay branching radios. Our results for 2 datum-fits are Br(p — ► 
7r + 7r~7r )| ( pp : v (7r o jK)) = (2.0 ± 1.64) x 10~ 2 and Br(p -> ir+ir- Tr a )\ {P p :V ^o , K ,r,,r,')) = (0-59 ± 
0.315) x 10~ 2 respectively, which are anomalously large and match each other within the errors. 

In order to pursue whether these anomalously large results of Br(p — > 7r + 7r~7r°) are rea- 
sonable or not, a theoretical estimation for p — uj— interference effects to the process of (p — > 
7r + 7r~7r°) has also been discussed in this paper. Following ref.QJ, we found that the contribu- 
tions due to so called hidden isospin-breaking effects are dominate for the process p — > 7r + 7r~7r°. 
The theoretical prediction is BR(p° — > 37r) ~ 0.2 x 10 -2 which is in good agreement with our 
datum-fit results. Then, considering this fact and noting that both result of r] — rf— angle 9 and 
the result of constituent quark ratio A = m u /m s obtained by the fits are also reasonable, we 
conclude that Br(p — > 7r + 7r~7r ) ~ 10~ 3 — 10~ 2 is reliable. 

Finally, we like to argue that in order to reduce the error-bar of Br(p — » ir + ir~ir°), more 
precisely experimental measurements to (J/ip — > PP, PV) are expected. The high quality data 
for J/ijj in the future BESIII would be useful. 

We would like to thank Zheng Zhi-Peng, Shen Xiao-Yan, Zhu Yucan, Yuan Chang-zhcng, 
Fang Shuang-shi for helpful discussions. This work is partially supported by NSF of China 
90103002 and the Grant of National Laboratory at the Institute of High Energy Physics, Beijing. 
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